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Abstract 

I introduce a new type of continuous-time quantum walk on graphs called the quantum 
snake walk, the basis states of which are fixed-length paths (snakes) in the underlying graph. 
First I analyze the quantum snake walk on the line, and I show that, even though most states 
stay localized throughout the evolution, there are specific states which most likely move on the 
line as wave packets with momentum inversely proportional to the length of the snake. Next I 
discuss how an algorithm based on the quantum snake walk might potentially be able to solve 
an extended version of the glued trees problem which asks to find a path connecting both roots 
of the glued trees graph. No efficient quantum algorithm solving this problem is known yet. 

1 Introduction 

In 1994, Shor's ground-breaking paper presented an efficient quantum algorithm for integer factoring 
[19]. Many other very important results have followed, but the true potential of quantum computers 
is still unknown. One step towards understanding it is to consider what are the differences between 
quantum and classical computers in a setting where we are given the input via a black-box oracle. 
Because we cannot know the inner workings of the oracle, it is easier to prove lower bounds in this 
model and therefore in many cases it is easier to show a separation between the best quantum and 
classical algorithms. 

Ever since Deutsch presented a simple oracle problem that can be solved on a quantum computer 
using fewer oracle queries than on any classical computer [10], scientists have tried to come up with 
more and more problems for which quantum computers outperform their classical counterparts; 
some very artificial, some quite natural. Bernstein and Vazirani gave the first example of an 
oracle problem which can be solved in polynomial time on a quantum computer, but requires 
superpolynomial time on a classical computer [3]. Shortly after that Simon gave an example 
in which this separation is exponential [20]. Yet it is unclear what are the best methods for the 
construction of efficient quantum algorithms. While quantum Fourier sampling is probably the most 
popular such method so far, many algorithms are also based on the concept of quantum walk. In 
particular, continuous-time quantum walks on graphs, introduced by Farhi and Gutmann [13], give 
rise to fast algorithms for NAND tree evaluation [11] and unstructured search [8,12]. Continuous- 
time quantum walks are also known to be able to perform universal quantum computation [6] , and 
they can solve some oracle problems exponentially faster than any classical algorithm [7,9]. 

In this paper we introduce a new type of continuous-time quantum walk on graphs, the basis 
states of which are not vertices of the graph, but paths in it of a fixed length n. We define this 
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new walk on an unweighted undirected graph G as a regular continuous-time quantum walk (the 
one introduced by Farhi and Gutmann) on a more complex weighted undirected graph Gn which 
is constructed from the original graph G. Since, to the author of this paper, the discrete-time 
classical counterpart of this walk in some sense resembles the behavior of a 'snake' of length n 
which is placed on a graph along its edges and which makes random decisions in which direction 
to go next, we call this new walk the continuous-time quantum snake walk. 

The main motivation for the continuous-time quantum snake walk is an extended version of the 
glued trees problem for which no efficient quantum algorithm yet is known. Childs et al. introduced 
the glued trees graph, consisting of two complete binary trees of the same height which are connected 
by a random cycle that alternates between the leaves of the trees [7]. Given a glued trees graph 
via a black-box oracle and the label of the root of one tree, the glued trees problem is to determine 
the label of the root of the other tree. While there is an efficient quantum algorithm based on the 
continuous-time quantum walk which solves this problem, no efficient quantum algorithm which 
finds a path connecting these two roots is known. We hope that the continuous-time quantum 
snake walk might lead to a quantum algorithm efficiently finding such a path. 

In Section 3 we analyze in detail the continuous-time quantum snake walk on the line. Even 
though most states stay localized throughout the quantum evolution, there are specific states which, 
under one reasonable assumption, move on the line as wave packets with momentum inversely 
proportional to the length of the snake. This motion is similar to the regular continuous-time 
quantum walk on the line. In Section 4 we discuss an algorithm based on the continuous-time 
quantum snake walk that solves the extended version of the glued trees problem mentioned above, 
that is, finds a path connecting the roots of the glued trees graph. However, it is not clear whether 
this algorithm runs in polynomial time. 

Independently from this work, a recent paper by McGettrick defines a discrete-time walk on 
the line similar to the continuous-time quantum snake walk on the line considered here [16]. 

2 Definitions and motivation 

In this section we define the main concept of this paper, the continuous-time quantum snake 
walk. The main reason for considering this new type of quantum walk is its potential algorithmic 
applications. Because of that we also define an extended version of the glued trees problem for 
which no efficient classical algorithm exists and no efficient quantum algorithm is known, but which 
might be solved efficiently using an algorithm based on the continuous-time quantum snake walk. 

2.1 Continuous-time quantum snake walk 

We will define the continuous-time quantum snake walk on one graph as the continuous-time 
quantum walk on another specific graph. Therefore, we define the continuous-time quantum walk 
first. Let K = (U, w) be a weighted undirected graph, where C/ is a set of vertices and w is a 
weight function, which assigns a weight w{ui,U2) G M to every pair (^1,^2) S U"^- Let C*^ be a 
Hilbert space having an orthonormal basis {\u) : u G U} called the standard basis, and let Hx be 
a linear operator acting on such that {ui\Hk\u2) = w{ui,U2) for all ui,U2 G U. Because K is 
undirected with all weights being real, Hk is Hermitian. The continuous-time quantum walk on K 
is defined as a quantum evolution in the space governed by the Hamiltonian Hk according to 
the Schrodinger equation. That is, if we fix an initial state of the walk to be |x(0)) ^ C^, then we 
can look on the walk as a function which maps time t > to the state = e~'^^*|x(0)) G C^. 

Now suppose G = {V, E) is an undirected unweighted graph, where y is a set of vertices and 
£^ C is a set of edges. Let Sn{G) be the set of all paths in G which have length n. Here we 
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assume that a path can visit a vertex multiple times. Let us call an element of Sn{G) a snake of 
length n. That is, a snake of length n is a vector s = {vq, . . . , Vn) G V^'^^ such that {vi-i,vi) G E 
for all / G [1 .. n]. Note that we assume all paths are directed {{vo,vi, . . . , Vn) and Un-i, ■ • • , vq) 
are not the same snake), even though the underlying graph G is undirected. 

Let s = {vq, . . . , Vn) G Sn{G) be a snake of length n. We say that s can move forward to a snake 
t G SniG) if there exists a vertex Vn+i G F such that t = {vi, . . . , Vn+i), and we write rnf{s, t) (we 
think of as a predicate). Similarly, we say that s can move backward to t G SniG) if there 
exists V-i G 1^ such that t = (w-i, . . . ,fn-i), and we write mb{s,t). We also consider mf{s,t) as a 
binary function on 5„(G)^ taking the value 1 if and only if the predicate mj is true for a pair (s, t), 
and similarly for ruh- Let as^t = fnf{s,t) +mb{s,t) G {0,1,2}, and let A„(G) be a matrix whose 
rows and columns are labeled by the elements of 5„(G) and An{G)s,t = as,t for all s,i G Sn{G). 
Clearly mf{s, t) if and only if mh{t, s), and therefore An{G) is a symmetric matrix. We can look on 
An{G) as the adjacency matrix of a graph having Sn{G) as the set of vertices, and possibly some 
edges having weight 2 instead of 1 (as a matter of fact, An{G)s,t = 2 if and only if there are two 
adjacent vertices u and v such that s = {u, v, u,v, . . .) and t = {v, n, v,u, . . .)). 

Definition 1. Let G„ be a weighted graph with the set of vertices S„(G) and the matrix of weights 
An{G). The continuous-time quantum snake walk on the graph G is defined as a continuous-time 
quantum walk on the weighted graph G„. 

Because we consider only continuous-time walks, we often refer to continuous-time quantum 
walks and continuous-time quantum snake walks, respectively, as regular quantum walks and quan- 
tum snake walks. Sometimes we use the word snake to refer, instead of an element of Sn{G), to 
the content of a quantum register corresponding to the Hilbert space C'^"^'^^ In some sense, we 
can think of the snake as being a path which can move form one position in the graph to another. 
In what sense the word snake is used should be clear from the context. 

2.2 Extended glued trees problem 

In order to provide an example of a quantum algorithm based on a quantum walk that can solve a 
certain black-box problem exponentially faster than any classical algorithm, Childs et al. introduced 
the glued trees graph [7] . The glued trees graph consists of two complete binary trees of height A^, 
which are connected by a random cycle that alternates between the leaves of the two trees. An 
example of the graph is shown in Figure 1 (for = 3). Let us also call N and the roots of both 
trees, respectively, the height and the roots of the glued trees graph. 

Suppose we are given a glued trees graph via a black-box oracle (given the label of a vertex 
as an input, the oracle outputs labels of all its neighbors) and the label of one root ri. The glued 
trees problem is to determine the label of the other root r2 (see Figure 1). Childs et al. show 
that, if we start a continuous-time quantum walk at the root ri, the walk quickly (in poly time) 
traverses the graph and ends up in a superposition state which has a large overlap on the root 
r2- We can implement continuous-time quantum walks on sparse graphs in the quantum circuit 
model efficiently (see [1,4,5]), therefore there is an efficient quantum query algorithm solving the 
glued trees problem. It is also shown in [7] that no classical query algorithm can solve this problem 
efficiently. 

Definition 2. Given a glued trees graph via a black-box oracle and the label of one of its two roots 
ri, the extended glued trees problem is to find a path in the graph connecting ri to the opposite 
root r2- 
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Figure 1: A glued trees graph of height 3. 



As far as I know, there is still no quantum algorithm known which solves the extended glued 
trees problem efficiently. It is not even known if such an efficient algorithm exists [5, Section 5.4]. 
However, the hope is that the quantum snake walk might lead to an efficient algorithm. 

3 Quantum snake walk on the line 

As an example of a continuous-time quantum snake walk, let us consider the quantum snake walk 
on the line, that is, the graph G = {V, E), where V = Z is the set of vertices and = {(x, x ± 1) : 
X G Z} is the set of edges. This example is relatively simple compared to quantum snake walks on 
other graphs. Nonetheless, understanding this walk later helps us to analyze snake walks on more 
complex graphs. 

3.1 The Hamiltonian 

It is convenient to think of the line also as the X-axis. For every snake {vq, . . . ,Vn) G V^^"^ on G 
of length n let X = vq be the start vertex of the snake and for I G [1 .. n] let ji G {0, 1} be such that 
vi = vi^i — {—ly'. This gives a one-to-one relation between the set of snakes 5'„(G) and the set 
Z X {0, 1}". From now on let us consider any snake s to be given as a pair (x, ji ■ ■ ■ jn) G Z x {0, 1}", 
and let |x)|ji . . . j„) denote \s). 

Every snake on the line can move forward to two other snakes and move backward to two other 
snakes. If a snake (x, ji . . . jn) moves forward, then the start vertex of the new snake is determined 
by ji. Its end vertex can either move in the positive or the negative direction of the X-axis. If a 
snake moves backward, the direction in which its start vertex moves is opposite to the direction of 
the first edge of the new snake. In other words, for every j G {0, 1}" we have 

mf{{x,ji . ..jn), {x - {-l)^\j2 . . . j„0)) and mf{{x,ji . ..jn), (x - {-l)^\j2 ■ • -Jnl)), 
rribiixji . ..jn), (x - 1, Iji . ..jn-i)) and m;,((x, ji . ..jn), (x + l,Oji . . . 

This allows us to obtain the Hamiltonian Hn = An{G) governing the quantum snake walk on the 
line. After some derivations we get 
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where 



and 



(1) 



ie{o,i}"-i 



+ E |0J')0'0| + |Oj)(ii| + Uo)(ij| + |ji)(ij| ) 

j6{0,l}"-i 

(see [18, Section 3.1] for details). For every fci, /c2 G M we have (/c2|fci) = ^(^i — ^2 mod 27r), where 
A is the Dirac delta function (see [6]). Hence, for any value A;, if jV') is an eigenvector of iin,k with 
the corresponding eigenvalue A, then \k) ® is an eigenvector of Hn corresponding to the same 
eigenvalue A. 

Let us focus on finding the eigenvalues and eigenvectors of -ffn,fc for an arbitrary h. This task 
becomes much easier if we express lin,k m a different basis. Consider two pairs of orthonormal 
vectors: 

ko,fc) = ^ (e-'» + e'^^ll)) and |ni,fc) = (e-'^^lO) - e^^^ll) 
bo) = ^(|0) + |l)) and |^i) = -L(|o)-|l)). 

For m G [0..2"' — 1] let us define a unit vector Im^) as follows. First, let |Ofc) = — il^o./c)*^"- For 
m G [1 .. 2" — 1], let m written in binary using [log2(rra)J + 1 bits be l?n'[iog2(m)J • • • "^i- Then for 
m G [1 .. 2" — 1] we define 

As an example,^for n = 3 we have |Ofc) = -i|Mo,fc)|'Wo,fc)K,fc), = |wo,fc)K,fc)l'Wi,fc), |2fc) = 
|^*o,fc)|wi,fc)|'yo), |3fc) = |^to,fc>|^ii,fc>|^^i>, |4fc) = |'Ui,/fc)|vo>bo), |5fc) = |ui,fc)bo)bi), |6fc) = |ui,fe)|vi)|'uo) 
and |7fc) = \u\^]yj\v\)\v\) . For any n G N and A; G M the set Bn^k = {|Ofc), . . . , |2" — 1^)} is an 
orthonormal basis of C''-'^'^-''", and [18, Section 3.1.1] shows that Hn,k expressed in this basis is 

/7„,fc = 4cosA:|0fc)(0fe|+2sinA:(|U)(0fe| + |0fe)(Tfc|) + 2 ^ {\2^k){mk\ + \mk){2^k\)- (2) 

m=l 

From (2) we see that {mk\HnM\lk) = whenever m and I have different greatest odd divisors, 
where we assume that the greatest odd divisor of is 1. Therefore we can block diagonalize i?n,fc 
with respect to -B^.fc into 2"~^ blocks. To be more precise, let us define 2"'"^ orthogonal projectors 

ni,fc = \dk){Ok\ + Ei=o \^'k){2h\ and Ui^k = E^^Zl'"^''^ \r2^k){l^^k\ for odd / G [3 .. 2" - 1]. We 
have Ill-^^^kHn,k^l2,k = whenever li / 12- Therefore, in order to get the full eigenspectrum of Hn^k, 
we consider the eigenvalues of each block separately. 
For odd Z G [3 .. 2" - 1] the block 

n-flogjil-l 

II/./,.//../ II'./, = 2 Yl {\f^h){i^h\ + \r2h){i^h\) , 

3=0 
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is basically twice the adjacency matrix of the line segment of length n — [log2 l] , whose eigenvalues 
and eigenvectors are well known (for example, see [2]). Therefore the eigenvalues of this block 
are fc-independent and the eigenvectors depend only on k due to the fact that the basis Bn^k is 
/c-dependent. However, because this dependence is relatively simple, as shown in [18, Section 3.1.3], 
the following lemma holds. 

Lemma 1. Let xi,X2 G Z and ji,j2 G {0,1}". If \xi — X2I > 2n, then {x2, j2\e~^^"^\xi, ji) = 
(a;2,j2|e"'^"*|xi, ji), where Kn = j^'^ \k){k\ ®Iii^kHn,k^i,k^k. 

The essence of Lemma 1 is that we can restrict our attention to the states in the space corre- 
sponding to the projector J^^ |^)(^| ® ^i,k dA; if we want to have an initial state which moves on 
the line further than distance 2n, which for algorithmic interests seems to be a reasonable request. 
The eigenvalues and eigenvectos of Hn restricted to this subspace, namely is determined by 
the eigenvalues and eigenvectors of the block Jli^kHn,k^i,k- 

Let {|1), . . . , |n + 1)} be a fixed set of n + 1 orthonormal, fe-independent vectors. Consider a 
linear fe-dependent isometry 

n 

Un,k = Y,\^k){y\ + |Ofc)(^^|. 

y=l 

We have Ili^kHn,k^i,k = Un,k^n,kU* f^, where 

n-l 

^n,k = 2^{\y +l){y\ + \y){y + l\) + 2sin/c (|n + l)(n| + |n)(n + 1|) + 4cos /c |n + l)(n + 1|, (3) 

y=l 

which is the adjacency matrix of the graph given in Figure 2. 

2 2 
• • • 

|T> |2) |3) 

Figure 2: The graph corresponding to the adjacency matrix ^n.k- 

If A: = mod vr, then sin k = and thus ^n,k is the sum of two orthogonal operators 
2 X^y=i (|y + + + and 4cosA;|n + l)(n + 1|. The eigenspecturm of the former op- 
erator is 2 times the eigenspectrum of the adjacency matrix of the line segment of length n — 1, 
and the eigenspecturm of the latter is trivial. Therefore, in the case when A; = mod vr we can 
get all the eigenvalues and eigenvectors of ^n,k- for every p G {^^j • • • > the unit vec- 
tor ^J'^^^Yly=l^^^yp\y) is an eigenvector of ^n,k with the corresponding eigenvalue 4cosp, and 
|n + 1) is an eigenvector corresponding to the eigenvalue 4cosA; (i.e., ±4). 

Unfortunately, in the general case we do not even know if there are closed form expressions for 
the eigenvalues of ^n,k- However, the following lemma, which is proved in Appendix A, gives us 
some useful information about the eigenvalues and eigenvectors of ^n,k for an arbitrary k. 

Lemma 2. Let us fix n gN and k ^0 mod vr. The equation 

2(cosp — cos A;) sin((n + l)p) = sin^ A:sinnp (4) 

has n + 1 distinct solutions in the interval (0,7r), and, if p is a solution of this equation, then 

X]y=i sinyp|y) + '^"^sinfc^^^ l^ ^) (unnormalized) eigenvector of ^n,k with the corresponding 

eigenvalue 4 cos p. 



rt-1 ) 



2 sin {k) 



n) |«+1) 




4 cos (k) 



6 



Corollary 3. The eigenvalues of ^n,k o'^e (^^l distinct for any value of k, and they all depend on 
the value of k. 

We refer to (4) as the p- equation. 



3.2 Even n and the median eigenvalue 

We do not know how to solve the p-equation, the solution of which would give us the full eigenspec- 
trum of Hn^k- Despite that, the p-equation allows us to obtain good approximations of eigenvalues 
and their derivatives. But let us first show that we can treat /c-dependent eigenvalues of Hn^k as 
infinitely differentiable functions of k. 

A complex-valued function is holomorphic if it is complex-differentiable in a neighborhood of 
every point in its domain. Holomorphic functions are known to be infinitely differentiable, and 
the class of holomorphic functions include polynomials, the exponential function, sine and cosine 
[15]. The notion of holomorphic functions can be generalized to vector- valued and operator- valued 
functions in an obvious way. From [14, Chapter II, §6.2] we have the following lemma 

Lemma 4. Let X he a finite complex Euclidean space of dimension m. Consider a holomorphic 
operator-valued function T which maps complex numbers to linear operators over X such that 
T(k) is Hermitian for all k £ M. Then for fc E M there exist a family of orthonormal basis 
{\(j)i{k)) : I £ [l..m]} of X consisting of eigenvectors ofT{k) and a family {Xi{k) : I £ [l..m]} 
consisting of eigenvalues ofT{k) such that Xi{k) and \(j)i{k)) are holomorphic functions of k for all 
I e[l..m]. 

Let A; (A;) be the l-th. largest eigenvalue of ^n,k- Because all the eigenvalues of ^n,k are distinct 
and ^n,k is holomorphic in k (which can be easily seen from (3)), according to Lemma 4, \i{k) 
is a holomorphic function in k for all / G [1 .. n + 1]. In order to give some intuition about the 
eigenspectrum of ^n,k, in Figure 3 we show how the eigenvalues of ^n,k and their derivatives 
depend on k G [0, 2tt] in the case when n = 8. These plots are obtained via numerical computation. 




Figure 3: The eigenvalues of ^n,k (on the left) and their derivatives (on the right) for n = 8. 

Figure 3 suggests that the derivative of any eigenvalue is if and only if = mod vr. The 
following theorem, which is proved in [18, Section 3.2.1] using the p-equation, confirms this obser- 
vation. 

Theorem 5. For every I G [1 ..n + 1] the eigenvalue function Xi{k) is strictly decreasing in the 
interval (0,7r) and strictly increasing in the interval (7r,27r). 
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For the rest of the section we consider only the case when n is even because it makes some 
calculations easier. Also, out of all n + 1 fc-dependent eigenvalues of Hn^k, let A(A:) be the median 
one, that is, the ^^^-th largest; we consider only this particular eigenvalue. The reasons for 
considering this particular eigenvalue are the fact that it can be well approximated and it seems to 
have the largest derivative, which most likely is good for a construction of states moving fast on 
the line. 



Note that the value of X{k) depends on n. Let A(A;) = 4arctan 



2 cos k 
sin*^ k 



and let us use a prime 



to denote derivatives with respect to k. In Appendix B we prove the following lemma. 



Lemma 6. For every n, X'{k) is bounded between ^^'^ {\ 



and 



A'(fc) 



1 + 



and we have 



A"{k) 



£ O (^) . Also, there exists no G N such that for all n > hq we have X"{k) = if 
^ mod vr. 



A"(fc) 

and only if k 

Corollary 7. The value of X{k) is bounded between (l 



) and^{l + l). 



As n increases, the function n\{k) and its derivatives nX'{k) and nX"(k) converge to the func- 
tions A(A;), A'{k) and A"{k), respectively. In Figure 4 we plot A(/c) and its first two derivatives to 
get some intuition about the behaviour of this 'eigenvalue' function. 




Figure 4: A{k) = 4arctan 
dotted lines, respectively). 



2 cos k 
s'm^ k 



(solid line) and its first and second derivatives (dashed and 



3.3 A localized superposition of 'median' eigenvectors evolving as wave packets 

Now that we have established some properties of the median eigenvalue X{k), let us consider a 
quantum snake walk on the line having an initial state |r/o) which is a uniform superposition of 
eigenvectors of corresponding to X{k). We will choose the phases of these eigenvectors in a 
particular manner so that |r/o) seems to be highly localized. Then we will show that, under one 
reasonable assumption, for n large enough the state |?7o) evolves as a left-moving and a right-moving 
wave packet, each propagating with the same constant momentum. But let us first give a precise 
definition of \t]q). 

Lemma 2 implies that there exists a unique vector- valued function !</>(•)) : M — )■ M"^^ such that 
^n,fc|0(fc)) = HkMik)), {(l){k)\(t){k)) = 1 and (T|(/>(A;)) > for ah E M. And Lemma 4 implies 
that \4>{k)) is continuous in A;, and therefore so is \'4'{k)) = Un,k\4>{^)) ■ Hence, \k) ® \'^{k)) is an 
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eigenvector of i/„ corresponding to the eigenvalue A(A;). For every x S Z define 

\V.) = ^ r e~"''\~k)(E>\m)dk. (5) 



Notice that \r]x) = S^ez Iv + ^){y\ ' \Vo)- The state |?7o) obviously belongs to the subspace V 
span( { (g) \^p{k)) : fc G [0, 2tt] } ) and remains in V under the evolution of Hn- Because {r]x\r]y) 



6x,y, where 6 is the Kronecker delta function, and 



\k)®\^{k))dk = \i)®\^{l)) 



for any / G M, {|?7x) : x G Z} is an orthonormal basis of V. 

Numerical results suggest that \r]x) is highly localized 'around the value x on the X-axis' (see 
Appendix C). For example, for n = 14, if we measure \rix) in the standard basis, with probability 
0.62 we obtain a snake starting at the position x it 1, with probability 0.26 a snake starting at 
X ± 3, with 0.09 a snake starting at x it 5, and so on (for even y we never obtain a snake starting at 
position X + y). Also, the probability of obtaining a snake starting at the position at least n units 
away from x seems to decrease exponentially with n, and the following hypothesis seems to hold. 

Hypothesis 8. For every y G Z and j G {0,1}" let Uyj = {y,j\rjx)- Then for every polynomial 
r{n) we have 

x—n +00 

Yl Yl \^v^i\ + Yl \^y'i\ ^ o{l/r{n)). 

Basically, Hypothesis 8 implies that for any \x) G V, if |(??a;|x)l is large, then for some y G 
[x — n + l..x + n — 1] the state | x) has large overlap with the space of snakes starting at the position 
y; meanwhile, if |(?7j/|x)l is small for all y G [x — n + 1 .. x + n — 1], then the state |x) has small 
overlap with the space of snakes starting at the position x. Therefore, the inner products between 
\x) and vectors of the basis {\rix) : x G Z} allow us to specify where on the line the state |x) is 
located. 

Let us analyze how the state |r/o) evolves under the Hamiltonian Hn in the limit of an asymp- 
totically large time of the evolution t. That is, for every a; G M let us consider (?/i^t|e~'^"*|r/o) for 
all t > such that uot G Z, and let us think of uj as the momentum. We have the following: 

{ri.t\e-'''-'\m) = ^ r e"("'^-^('=)) dk. (6) 

Let fu]{k) = uk — \{k). When t is large, e'*(^'^~'^('^)) = e^^^'^^^'^ rapidly changes as fuj{k) changes. 
This means that in an interval where fu){k) changes the contribution from adjacent subintervals 
to the integral (6) nearly cancels out [17]. Therefore the most contribution to the integral comes 
from values of k where the function fu){k) is stationary, that is, where /^(A;) = 0. The method of 
stationary phase approximation makes this statement more rigorous, and we will use it to estimate 
the integral (6) in the limit of asymptotically large t. 

From the p-equation one can get that A'(|) = —iSji ™d -^'(t") ~ nT2' -^^'^Hf ) = and 



' n+2 '^^^^ \ 2 I n+2 ' V 2 (n+2)=* 

A('^)(^) = — ^|^'^2)^3 ^ ■ Let us consider n > hq, where uq is the same as in Lemma 6. Lemma 6 implies 

that the derivative A'(A;) is strictly increasing in the interval (|, ^). For each u G (— 

let kcj be the unique point in the interval (|, ^) such that fui{kuj) = 0. Notice that for uj such 
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that > the function fio{k) has no stationary points. The method of stationary phase 
approximation (see [2,17]) applied to the integral (6) gives us that, in the limit of asymptotically 
large t. 



for every polynomial r{t); 
if w = then 



, -\H u \ -^t (n + 2)r(l/3) 1 



- n + 2 



27r -^4^/3(3^2 + 4) 



if'^e (-42'42)'then 

(r?.t|e-'^"*|r?o) « cos (t (fc^ - |) - A(A;^)) - |) 



As we can see, if we measure e~'^"*|7?o) in the basis {\r]u)t) '■ G Z}, most likely we will obtain a 
state |r/^t) such that A"(/c^) is close to 0. According to Lemma 6, A"(A;^) is close to when k^^ is 
close to I and ^ (see Figure 4) and therefore w, the momentum, is close to i^^- Hence, if indeed 
the state \r]i^t) is localized around the position tot on the X-axis, as Hypothesis 8 asserts, then we 
can see that for n > uq the state \t]o) evolves as a left-moving and a right-moving wave packet each 
propagating with momentum In Appendix C we present numerical data that suggest that 

this type of evolution can be already observed for small values of the time t and length n. 



4 Quantum snake walk on the glued trees graph 

In this section we discuss how the continuous-time quantum snake walk might potentially lead to 
an efficient algorithm solving the extended glued trees problem, which we defined in Section 2.2. 
Our aim is to come up with an algorithm which makes O(polyA^) oracle queries, where is the 
height of the glued trees graph and, without loss of generality, we assume each vertex is encoded 
using 9(A^) bit label. 



4.1 The main idea 

Let us consider the following expansion of the glued trees graph. Suppose G is a glued trees graph 
of height N, and ri and r2 are its roots. For M > N we construct an expanded glued trees graph 
G^^ as follows: take 2^^~^ instances of graph G, create two complete binary trees Ti and T2 each 
of height M — N, and then make bijective associations between leaves of Ti and instances of ri and 
between leaves of T2 and instances of r2 to obtain G^^ . An example of the graph for a given 
glued trees graph G of height 2 is given in Figure 5. Basically, G^ looks exactly like a glued trees 
graph of height M, but instead of one cycle gluing the trees, now we have 2^~^ cycles. 

For the rest of the paper we will focus on the continuous-time quantum snake walk on G*^. 
Given an oracle for G, since we know an efficient algorithm finding the label of r2 [7], we can 
efficiently simulate the quantum snake walk on G^^ for M G 0(poly A^) and the length of the snake 
n E 0(poly A^) (see [18, Section 4.1.3]). The sketch of an algorithm potentially solving the extended 
glued trees problem would be as follows: 



10 



1. choose M,n G O(polyiV) such that n > 2iV + 1; 

2. construct a particular initial superposition for a quantum snake walk on G^^ which contains 
only snakes on Ti (construction of any such state involves no oracle queries at all); 

3. run a quantum snake walk on G*^ for a certain amount of time t E 0(poly A^) and measure, 
thus obtaining a snake of length n; 

4. if the obtained snake contains two vertices such that one belongs to Ti and other to T2, then 
the snake has to go completely through one instance of G, and from this snake we can extract 
a path connecting ri and r2; otherwise return to Step 2. 




Figure 5: An expanded glued trees graph G for N = 2 and M = 4. One instance of the original 
glued trees graph G is highlighted using thick edges. 

First let us show that the analysis of the quantum snake walk on this relatively complex graph 
becomes much easier if we restrict the choice of initial state to certain superpositions. For x E 
[— M .. M + 1] let Yx be the set of all vertices in G^^ which are at distance M + x from the root of 
Ti (see Figure 5). Now, for x E [— M .. M + 1] and j E {0, 1}" let S{x,j) be the set of all snakes 
s = {vq, . . . ,Vn) such that vq E Yx and for all I E if E Yz, then vi E For 

S{x,j) 7^ let \x,j) be the uniform superposition over the elements of S{x,j), and let us refer 
to span{|a;,j) : S{x,j) 7^ 0} as the column subspace. Due to symmetry, the column subspace 
is invariant under An{G^), the Hamiltonian governing the quantum snake walk on G*^, and the 
action of An{G^^) on the column subspace does not depend on which cycle in particular glues the 
two binary trees in G. Also, let us call the subgraph of G^^ having the set of vertices Yq U Yi and 
all the edges originally connecting Yq and Yi the glued part. 

From now on let us consider M such that M ^ n. For |x| <C M and j E {0, 1}", if we run 
the quantum snake walk on G*^ starting from the state \x,j) for time small enough, then we will 
never get close to the roots of G^'^. Thus, in such a case the walk behaves almost exactly the 
same way on all G^ having M' > M, and for the sake of simpler analysis we assume M = 00 
(similarly as the regular quantum walk on a finite line segment, assuming the initial state is chosen 
to be far from the endpoints of the segment, can be approximated by the walk on the infinite line). 
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Therefore the Hamiltonian governing the quantum snake walk on G restricted to the column 
subspace effectively is 

(- max{|j|±,0) 
(2 ( k + 1, Oi) + \x, jl){x + 1, Oil ) 

+ ( |x - 1, J0| + |x,iO)(x - 1, lj\ ) (7) 

+ ^/2(|x + l,0j)(x,j0| + |2;,jO)(x + l,Oi|) 
+ \/2(|x- + l,li|)) 

-min{|j|±,0) 

+ (V2(|x + l,0j)(x,jl| + |x,il)(x + l,0j|) 

a;=l— max(|jr|±,0) 

+ V2{\x-l, lj){x,3^\ + |x, jO)(x - 1, lj| ) (8) 
+ If|j|±>o[2; 1] ( |x + 1, Oi)(x,iO| + |x,iO)(x + l,Oj| ) 

+ If|,-|±>o[l; 2] ( |x - 1, jl| + |x,il)(x - 1, lj| )) 

+00 

+ ((|x + l,Oi)(x,jl| + |x,jl)(x + l,Oj|) 

x=l-min(|j|±,0) 

+ 2 ( |x - 1, lj)(x, jO| + |x, jO)(x - 1, lj| ) (9) 
+ V2{\x + 1, Oj)(x, jO| + |x, jO)(x + 1, Oj| ) 

+ \/2 ( |x - 1, lj)(x, jl| + |x, jl)(x - 1, lil 

where |i|± = - ^^"^^(-1)^. and If,,-, ^>o[ai;a2] is ai if \j\± > 0, and 02 otherwise (see [18, Section 
4.2] for details). 

4.2 Quantum snake walk on the infinite binary tree 

Here let us consider the continuous-time quantum snake walk on restricted to superpositions 
\x,j) such that —M <C x <C —n. We care only about the first sum (7) in the expression for Hn, 
and, because x <C —n , similarly as before, we can approximate an action of Hn by one of 

f-2TT 

Hn = / \~k)Ck\'g)Hn^kdk, 
Jo 

where 

Hn,k = Y ( ^UO)(Oil + 2Ul)(0j| + |li)(jO| + V2\lj){jl\ ) 

i6{0,l}"-i 

+ Yl e-i'(^^|0j)(i0| + 2|0j)(il| + U0)(li| + ^/2|il)(li|). 

ie{o,i}— 1 

The analysis of Hamiltonian Hn and the quantum evolution governed by it will highly resemble 
the analysis of Hn above. Let us redefine \uo^k) = :^(\/2e~''^|0) + e''^|l)), \ui^k) = ~ 

V2e^^|l)), \vo) = 4^(|0) + V2\l)) and |t;i) = -j^{V2\0) - |1)), and let us define the orthonormal 
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basis {|Ofc), . . . , 12*^ — 1^)} as before except that now we are using these new definitions of |uo,A;) 
Ki,fc)) l^^o) ^-iid l^i) to define {ink) for each m £ [0 .. 2" — 1]. In this new basis 



Hn,k = 4\/2 cos A; |0y!c)(0fc| + (3sin/c + icosA;) |lfc)(Ofc| + (SsinA; — icosfc) |Ofc)(lfc 



+ 3 ^ {\2mk){mk\ + \fhk){2mk\). 



m=l 



Again we can see that only n + 1 eigenvalues of Hn,k are fc-dependent, and those are the eigenvalues 

of operator ^n,k = U* kHn,kUn,k, where lfn,k = ^y=i |2"~^fc)(l7| + |Ofe)(n + 1|. Along similar lines 
as for Lemma 2 we can prove the following lemma. 

Lemma 9. Let us fix n €zN and A; € M. The equation 



has n + 1 distinct solutions in the interval (0, vr), and, if p is a solution of this equation, then 6cosp 
is an eigenvalue of ^n,k- 



Similarly as for the quantum snake walk on the line, numerical results suggest that \ fix) is localized 
on the tree Ti at the distance |x| away from the glued part. If this property indeed holds, we can 
show that for n > ho and for asymptotically large time of the evolution t the state \fjx) evolves as 
two wave packets moving with momentum in the opposite directions, that is, one moving up 
the tree and the other down the tree. However, now we need to be careful when we talk about 
asymptotically large t, because we still want t to be small enough so that the wave packets do not 
get close to the root of Ti or the glued part. 

For the purpose of our algorithm, we are interested in only one of those two wave packets, the 
one which moves towards the glued part. Our hope is that in 0(poly A^) time this wave packet will 
propagate far enough to reach T2 and with high probability will give us a snake containing vertices 
from both Ti and T2. The regular quantum walk on the line started from a single vertex evolves 
as two wave packets moving in the opposite directions with momentum 2 [5, Section 3.3.2], but for 
any u € [0,2] we can create a wave packet which moves with momentum u in one single direction 
(see [6]). We would like to do an analogous thing for the quantum snake walk. For ko S (vr,27r) 



6(3 cosp — cos k) sin((n + l)p) 



(1+8 sin^ k) sin np 





Let \(l){k)) G C"+^ be the unique eigenvector of ^n,k corresponding to the eigenvalue X{k) such 
that {4>{k)\4>{k)) = 1 and (T|(^(/c)) > 0, let \ij{k)) = C/„,fc|0(A;)), and let 
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and fj > consider the state 

\Cxo,ko,cr) - 




-ikxQ 



\k) ® \ij{k))dk 



Z&7, 



eW7^ s)fj (erf + iz(7)) 



\'nxQ+z 



1 

E 



1 

2^ 



>2L^J 



where erf(-) denotes the error function and K(-) denotes the real part function. One can show that 
(?xo,fco,cr|'?xo,fco,o-) = 1- Numeric results suggest that |i^a;o,fco,o-) evolves as a Gaussian-shaped wave 
packet moving towards the glued part of G^^ with momentum A'(/co) (see Appendix D). 

For any snake in the superposition e"^^"*!^^;^^^^^^), other than the position of its initial vertex, 
we are also interested in how many different depths in the tree T\ its vertices possesses. To be more 



precise, for a snake s = (uq, 



let (zo, . . . , z„) G Z"+^ be such that vi € Y;, for alU G [0 .. n]. 



We define the span of s to be (min; z;, max; zi) and the span length of s to be max; zi — min; zi. We 
are interested in snakes with reasonably large span length, because, after all, if we obtain a snake 
which contains a solution of the extended glued trees problem, the span length of this snake will 
be at least 2N + 1. Numerics presented in Appendix D suggest that, if we measure e~^^"^\^xo,ko,(T) 
in the standard basis and output the span length of the obtained snake, then the expected value 
of this random variable is in i}{^/n), which is good. However, the wave packet might drastically 
loose its average span length when it gets close to the glued part and the action of the operator Hn 
cannot be approximated anymore by one of Hn- In Appendix E we describe numerical data which, 
unfortunately, suggest that it indeed may be the case. 



4.3 Scattering on the glued part 

Let us consider yet another analogy to the regular quantum walk on lines. Consider the regular 
quantum walk on a graph A constructed from some finite graph B by attaching two semi-infinite 
lines to two vertices of B. On either of those two lines we can construct a wave packet which moves 
towards the graph B. What are the probabilities of this wave packet being reflected form B and 
it being transmitted through B (to the other line) can be calculated by inspecting the eigenvalues 
and eigenvectors of the adjacency matrix of A and then using the standard scattering theory (see 
[6], for example). 

Similarly, for the quantum snake walk on G^'^ we can construct wave packets which move on Ti 
towards the glued part (to be more precise, we seem to be able to construct such wave packets, since 
we have not rigorously proven that yet), we can also construct wave packets which move on Ti away 
from the glued part, and we can do the same on T2. In order to understand what happens to the 
wave packets moving towards the glued part once they reach it, let us inspect particular eigenvectors 
and eigenvalues of Let \ilj{k)) be the unique vector satisfying (j|^(A:)) = {^l>{k)\j © 1") for all 
j E {0, 1}". Numerical data for n up to 10 agree with the following hypothesis. 

Hypothesis 11. For every A; € M A(A;) is the eigenvalue of Hn coresponding to the eigenvector 

— n+l n— 1 +00 

^ |x) © (e'^1V^(A;)) + R{k) e-'^nv^(-fc))) + J] |x) © a,,,(fe)|j> + ® T{k) e'^^im), 

x=—oo x=—n+2 jG{0,l}" x=n 
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where axj{k) G C, 



1 - 2e2i*^ ^ _ V2(e-2ifc _ 3 + 2e^''') 



Let us denote the eigenvector given in Hypothesis 11 by \fi{k)). Assuming Hypothesis 11, the 
standard scattering theory suggests that the wave packet (started as) ICxo.fco.o-)) where ko £ {it,2-k) 
and a is assumed to be small, will be reflected from the glued part with probability \R{ko)\'^ = 
^^sL^ko ^^"^ transmitted through it with the probability |r(A;o)|2 = ^g^- Also, arg'(r(A;o)), 
called the effective length, divided my the momentum of the wave packet describes for how much 
time the glued part delays the wave packet propagating through it. It is really interesting that 
the transmission coefficient T{k) does not depend on n, the length of the snake, and it still stays 
the same if we consider eigenvectors of Hn corresponding to any other /c-dependent eigenvalue of 
Hn^k, not only the median eigenvalue. Even more, already in [7] Childs et al. proved that for the 
regular quantum walk on the probability of a wave packet corresponding to the vector | k) being 
transmitted through the glued part is yJI^- Figure 6 shows that the wave packet \(,xo,ko,cT) has 
the highest probability of being transmitted through the glued part and has the minimal effective 
length of the glued part for fco = It also seems that A'(/co) reaches its maximum at /cq = ^ 

(we have A'(^) = |^), therefore the faster a wave packet is, the most likely it will propagate 
through the glued part (note: this is the case only for wave packets corresponding to the median 
eigenvalue) . 




Figure 6: The likely probability of the wave packet l^^xo.fco.o-) propagating through the glued part 
on the left, and, on the right, its corresponding effective length of the glued part. 

Even if we rigorously proved that we can construct localized wave packets having momenta as 
stated above, the span lengths of those wave packets are as stated, and so are the transmission 
probabilities, such a proof most likely would use the method of stationary phase approximation, 
which deals with the case when all (i.e., initial, reflected and transmitted) wave packets are far 
away from the glued part. But the problem is: what happens during the time when wave packets 
encounter the glued part is exactly what we are interested in for the analysis of our algorithm. 
Unfortunately we cannot obtain enough numerical data to have a good intuition what is the maximal 
probability obtaining a snake connecting Ti and T2 for any given N and n > 2N + 1. Also, n 
necessary for an efficient algorithm may be very large because we only restrict it to be in 0(poly N). 
Nonetheless, a numerical inspection of the vector |^(fc)) done in Appendix E suggests that the 
expected span length of snakes in a wave packet is significantly reduced when the wave packet 
propagates through the glued part. We currently do not have any candidates for a pattern this 
reduction follows, so it is not clear if the reduction is exponential or not. 
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5 Discussion 



In this paper we introduced the continuous-time quantum snake walk and discussed one of its 
potential algorithmic applications, the extended glued trees problem. Assuming the algorithm 
sketched in Section 4.1 which uses wave packets as initial states indeed efficiently finds a path 
connecting both roots of the glued trees graph, it almost certainly will not find a shortest path. 
It would be interesting to see if there exists an efficient algorithm finding a shortest path, or an 
efficient quantum algorithm for any other pathfinding problem, which is based on the quantum 
snake walk. 

On the other hand, if the algorithm we propose is not efficient, despite that the wave packets 
(most likely) are being transmitted through the glued part with high probability and in polynomial 
time, it would suggest that possibly there are some fundamental limitations which restrict us from 
solving the extended glued trees problem efficiently. Finding and understanding such limitations 
would be an interesting and useful thing to do. 

Another potential direction for future research is to consider discrete-time analogues of the 
continuous-time quantum snake walk (for example, by generalizing McGettrick's results [16]) and 
their applications. 
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A Proof of Lemma 2 

Let \(p) = J2y=i '^ylV) be an eigenvector of ^n,k and A be its corresponding eigenvalue. The 
maximum absolute column sum norm of (1) is at most 4 (it is exactly 4 unless n = 1), which 
implies that all the eigenvalues of Hn,k are at most 4 by the absolute value, and therefore so are 
the eigenvalues of ^n,k- In fact, it is easy to see that |A| = 4 would imply k = mod it and 

oc |n + 1), which is also the only case when oi = 0. Hence, there exist unique p £ (0,7r) and 
c 7^ such that A = 4cosp and oi = csinp. 

We have 

Aai = A(T|0) = (T|$„,fe|(/>) = 2{2\ct)) = 2a2, 

from which we get 02 = (4cosp • csinp)/2 = csin2p. Now let us use induction. Let 2 < I < n, and 
let us assume Uy = csinyp for all y G [1 .. /]. We have 

Xai = Xm = = (2(1^1 + 2{TTim = 2a,_i + 2ai+^, 

which gives us 

ai-^-i = (4cosp • csinlp — 2csin(/ — l)p)/2 = csin(/ + l)p. 
Therefore, by induction, Uy = csinyp for y £ [1 ..n]. Similarly, 

Aa„ = A(n|(/)) = {n\^n,k\<t>) = (2(n - 1| + 2sinfc (n + l\)\(f>) = 2a„_i + 2a„+i sink 
implies On+i = c ^'^^sink^^ ■ Finally, 

Aa„+i = A(n + = (n + l|<I>„^jt|(/)) = (2sinA:(re| + 4cosA;(n + l\)\cj)) = 2sinfc • a„ + 4cosA; • a„+i 
gives us 

2(cosp — cos A;) sin((n + l)p) = sin^ ksinnp. (10) 

One can see that (10) is both the necessary and sufficient condition for 4 cos p to be an eigenvalue 
of Also, as shown above, an eigenvalue of ^n,k uniquely (up to a global factor) determines 

the eigenvector corresponding to it. This implies that all the eigenvalues of ^n,k are distinct, thus 
(10) has n + 1 distinct solutions in the interval (0,7r). 
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B Proof of Lemma 6 



For n = 2 we have A(A:) = 2 cos(A;) and the lemma holds, thus let us assume n > 4. We know 
that A(0) = 4 cos and A(7r) = 4 cos ^^^-^ , therefore a corollary of Theorem 5 is that |A(A;)| < 
4 sin 2n+2 • P(^) ~ f ~ ^ \- 2n+2 ^ ^2n"+2^ ] such that X{k) = Acosp{k), and therefore 
e{k) G [-f , f]. Hence X{k) = 4sin(gi), where 6l(fc) is the solution of 

2 ^cos k - sin (^^^^^'('fc)^ ^ cos 6l(A;) = sin^ k sin , 
which is just a rewritten p-equation. By simple trigonometric derivations this gives us 

2 ^°^'-.^:l^) +sinfg§) 



cos 



\n+l J 



Notice that {tan9{k))' = 9'{k){l +tan^ ^{k)), and we already have an expression for tan0(A;). Thus 
we can obtain an expression for 9'{k), from which 



A'(fc) 



40'(fe)cos^ 16sinfe-cos2^ 

^ + 1 n(3 + cos 2k) +4- 4(n + 1) cos /fc • sin ^ ' 



Because A(k) = A(— k) and A(| + k) = — A(| — k) for any k G M, without the loss of generality let us 
assume k G (0, f ); hence O^k) G (0, f ). Therefore < sin g| < 1 - < cos < 1 

and 

< 4(n + 1) cos k ■ sin -^^^ < 2tt cos k < 2tt, 
n + 1 

which implies 

X!ik) > 16sinfc ^ 16 sin /c / 2 

' ^ - ~n(3 + cos2A;) - 2(7r - 2) " ~n(3 + cos2/c) V ^ 

^ < ISsinfe- (l-8(^) ^ 16sinfc / _ 2 
^ - n(3 + cos2fc) + 4 " n{3 + cos2k)\ n 

Since we already have expressions for X'{k) and 0'{k), we can get that 

9{k) I n(cos A; + cos^/c + 2 cos sin^fc) 



A"(A;) = 8cos^ 



^ + V (n{l + cos^A;) + 2 - 2(n + 1) cos k • sin 



2 



2(n - 1 + 2 cos^/c) sin ^ - 2 cos k 

7 2 (11) 

f n(l + cos2A;) + 2 - 2(n + 1) cos k ■ sin ^ j 

4(n + 1) cos k sin^k ■ cos'^ \ 



+ 



n(l + cos^A;) + 2 — 2(n + 1) cos • sin ^fpj j 



The first term of this equality is in while the other two are in O (^)- More precise analysis 
of the first term reveals that there is a constant a > such that 



X"{k) - ^ 



^ a 
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Regarding the second part of the lemma, along the similar lines as for \"{k), one can show that 



there also exists a constant 6 > such that 



A(3)(A;) 



A(^)(A:) 



< 



Because of that, there exist 



no G N and K e (0, f ) such that A"{k) + ^ < for aU /c G [0,k] and A^^\k) - ^ > for all 



b_ 
no 



k G [k, f]. Consider n > uq. We have X"{k) < for all k G [0,k]. The fact that A"(f) = and 



A(3)(fc) > for aU k G 
/cG [0,f). 



[k, f ] implies X"ik) < for A; G [k, f ). Thus, A"(A;) / for all n > no and 



C Numerics for the quantum snake walk on the hne 

Fix even n G Z. For any given k we can numerically calculate the median eigenvalue A(A:) of ^n,k 
and its corresponding eigenvector \(p{k)). Then we can get a numerical approximation of |?7o) by 
approximating the integral (5) defining it by a Riemann sum. 

For every unit vector \x) G C^^i"'!}", y G Z and j G {0,1}" let Xy,j = (yijlx) and let 
\x)x = T.j(z{o^iynXx,j\x,j). Then, in the standard basis, we have ||x)a;|i = I]je{o,i}" l^^.il ^"^^ 
llx)a;|2 = X^je{o 1}" IXi'jP- If we measure \x) in the standard basis, llx)^!! is probability of 
obtaining a snake which starts at the position x; therefore the probability distribution {||x)a;|2 • 
X £ Z} gives us good intuition about the 'location' of the state \x) on the line. 

Lemma 12. For every even x €z Z we have \r]o)x = 0. 

Proof. Let x be even, so that e''^^ = ei('=+'^)^'. From the definitions of |?7o) and \k) we can see that 
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Figure 7: For n = 14, the probability ||??o)a;|2 of obtaining a snake starting at the position x by 
measuring jryo) in the standard basis (the linearly scaled plot on the left and the logarithmically 
scaled plot on the right). Note: numeric values below 10"^*^ cannot be trusted because they 
approach the precision of the floating point arithmetic used. 
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As an evidence that \i]q) is locaHzed around the position on the Hue, in Figure 7 we present 
numericahy computed values of the probability ||??o)a;|2 the case when n = 14 for x close to 0. 
Lemma 12 states that this probability is for even x, therefore let us now focus on odd x. We can 
see that ||?/o)x|2 seem to decrease exponentially as |x| increases, and especially large jump seem to 
be from ||r?o)i4-i|i » 1-61 • 10"^ to \\vo)u+i\l » 1.47- lO^^. 

A similar pattern, a large gap between log(||?7o)n-il2) ^iid logd I%)n+il2)) was observed for 
values of n other than 14. This is the main motivation why Hypothesis 8 is stated in the way it 
is stated. The reason why we choose 1-norm in it instead of 2-norm is a technicality: it allows us 
to make stronger statements about superpositions of vectors \rix), where x ^ Z. Numerical data 
showed in Figure 8 support Hypothesis 8. 
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Figure 8: In the logarithmic scale, 1-norm of the part of the state |r/o) corresponding to snakes 
starting at least n positions away from (note: I770) is n-dependent itself), that is, J2x=-oo \\'no)x\i + 
^x=n \\'no)x\i- For n = 2 this value can be shown to be 0. 

Again, let n = 14, and let us compute how |7/o) behaves under evolution of Hu. Section 3.3 
states that, under certain assumptions, for asymptotically large time t the state \r]o) evolves as 
two wave packets each moving with momentum = ^. We can see in Figure 9 that this kind 
of motion can be observed already for smaller values of t. Notice that peeks of the probability 
distribution around positions ib|t get relatively higher for larger t. Higher probabilities around the 
origin x = are likely there due to the fact that , has a local maximum when X'(k) = 0. 

D Numerics for the quantum snake walk on the glued trees graph 

The same way as in the case of the quantum snake walk on the line, we can express er^^'^^\^XQ,kQ,fj) as 
an integral and then approximate it by a Riemann sum. Also, for every unit vector |x) G C^^'^^'^J'" 
^(j--Zx{o,i}" bging the column subspace) and x G Z let ||x)x|2 = (xl ( ® ^{0,1}") Ix); which, if 
we measure |x) in the standard basis, is the probability of obtaining a snake whose initial vertex is 
in Yx- Figure 10 shows that the state |^xo,fco,o-) evolves as a Gaussian-shaped wave packet moving 
with momentum X'{ko) towards the glued part. The diffusion of the wave packet can be explained 
by the Heisenberg uncertainty principle: we cannot know both the position and the momentum of 
the walk to an arbitrary precision. 

For every x E Z and j € {0, 1}" the span length of all snakes in S{x,j) is the same, and it 
depends only on j. Therefore we can easily construct the unique observable Qn acting on C^^'^^" 
such that for every \x) G C'^^^'^'^^" , if we measure \x) in the standard basis and return the span 
length of the resulting snake, the expected value of this random variable is <Xi Qn)\x)- For 
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Figure 9: The initial state |77o) evolved for time t = 400 and t = 800. 




Figure 10: The propagation of the wave packets \$,xo,ko,(T) for fco = x (*^® ^^^^ plot) and ko = ^ 
(the right plot) in the case when n = 8 and a = 1/20. The initial position of the wave packet is 
given with the dashed lines, and the position after time t = 100 with the solid lines. Numerically 
we get A'(^) 0.42 and A'(^) «i 1.13. 
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e '^"*|Cxo,fco,T) we have 

ert(7r/v2cr) Jfco-7r crV27r 
>mm(V;(A:)|Q„|V'(A;)). 

Figure 11 shows a plot of numerically obtained values of {ip{k)\Qn\ip{k)) for all values of even n 
up to 20. Prom this plot alone it is not completely clear how fast minfc(^(A;)|(5n|V'(^)) grows with 
n, but it seems to grow with the same rate as maxk{ilj{k)\Qn\4^{k)) . For any j E {0, 1}*^ we have 
{j\Qn\j) ^ |2| jIh — ?^|, where \-\h is the Hamming norm, therefore for j uniformly chosen at random 
we have {j\Qn\j) £ ^{\/n). Because 

i^( 2 )) = -^y^ E(i^i) ^ i^i))""' ® (i^o) ® I v)^"/'-', 

which is a quite balanced superposition over the vectors of : j G {0, 1}"}, we conjecture that 
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Figure 11: The span length of the snake |x) (8> \ipik)) for an arbitrary x G Z. 




E Evidence against efficiency of the algorithm 

Consider the walk on the expanded glued trees graph started as the wave packet \(,xo,ko,(T) , which 
moves towards the glued trees, and consider the maximal probability that, if at any point during 
the walk we make a measurement, we obtain a snake which contains vertices from both Y^m and 
Ym+i- We are interested to choose n £ 0{po\y N) such that this probability is at least inverse- 
polynomially large in N for m = N . It is not clear whether such choice is possible, and here 
we present numerical data suggesting that n we need to choose may me very large, if not even 
exponentially large in A^. 

We introduced the vector in Hypothesis 11 as an eigenvector of Hn- Unless its coefficients 

(y-x,j{k) behave very unexpectedly and ai^ {ax,j{k)) <^ —n, for <^ —n and a £ Q. the 

state \S,xo ^kQ,a) can be weh approximated by a superposition of states |/u(A;)), specifically, 

\U,k„a) cc / e-^^e-'''''^"\k)0\^P{k))dk^-= / e-^^e-"=-'°|/i(fc)>dA;. 
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Because is an eigenvector of Hn for every € M, the state e~^^"^\^xo,ko,<^) '^^^ ^^^^ 

approximated by a superposition of states too. Therefore, if for ah A; G M the state 

has a small overlap with the space of snakes connecting vertices from both YIat and YAf+ij so does 
e~'^^"-^\(^xo,ko,(7) ■ Now, let us numerically calculate the coefficients axj{k) and see what we obtain 
if we "measure" in the standard basis (note that has an infinite Euclidean norm). To 

be precise, we will be interested in the "probabilities" Px,aik) defined as follows. For even a < n 
and X E Z and for odd a < n and x G Z + ^ let Px,a{k) be the projector to the space of all snakes 
having span (x — + |) G Z^. Let px,a{k) = {/i{k)\Px^aik)\fi{k)) . We numerically calculated 
these values for n = 10 and Figure 12 shows for every a G [2 ..n] how Px,a{^) depends on the value 
of X (in this particular case px^i 

(^) = for all x). 

P.-,„(3jr/2) p,.aOJT/2) 
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Figure 12: The "probability" of obtaining a snake having span (x — |, x + |) if we "measure" 
\fi{'^)). Note: the scaling of the plots differs. 

A similar pattern as in Figure 12 is observed if we consider values of Px,a{k) for k other than 
^ or "measure" eigenvectors of Hn corresponding to other ^-dependent eigenvalues of Hn,k, not 
only the median eigenvalue. This means that the overlap between these eigenvectors and the space 
of snakes having a small span length drastically increases close to the glued part, and decreases 
for snakes having a large span length. Therefore, the expected span length of wave packets will 
get smaller when they come close to the glued part, and will increase again once they have been 
reflected from or gone through it. However, it is not clear what pattern this relative reduction of 
the span length follows when n increases. 
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